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^ Abstract. The interface between an unstable state and a stable state usually develops a single 
Q\ confined front travelling with constant velocity into the unstable state. Recently, the splitting of 
Q^j such an interface into two fronts propagating with different velocities was observed numerically 

in a magnetic system. The intermediate state is unstable and grows linearly in time. We 
q first establish rigorously the existence of this phenomenon, called "dual front," for a class of 
| 7 j structurally unstable one-component models. Then we use this insight to explain dual fronts for 

a generic two-component reaction-diffusion system, and for the magnetic system. 

<N 

t> 
o 
o 

m 
o 

+-» 

-I— > 

c3 
Oh 



-2- 

1. Introduction 

Parabolic differential equations, such as the (complex) Ginzburg-Landau equation [CE, CH], 
occur in many domains of physics. We are here interested in this equation when considered 
on an infinite domain. Then, a well-known phenomenon is the formation of fronts (travelling 
waves), which are solutions whose shape is fixed in a frame moving with speed c. For example, 
for the equation 

u = u" + u — u 3 , (1-1) 

one has the explicit solution 

u(x,t) = ^ (l -tanh(2" 3/2 (a; - 3 • 2~ 1/2 *))) , 

-I ley 

which moves with speed c = 3 • 2 1 . In a recent paper, a new phenomenon was observed 
[EBS]: A model for the dynamics of an easy-plane ferromagnetic material was considered, and 
it led to the surprising discovery that several concurrent front speeds can be observed, with 
a leading front connecting an unstable state to an intermediate unstable state, which in turn 
gives rise to a transition to a stable final state. As time goes on, the spatial region in which 
the intermediate state is present is seen to grow. We shall call this phenomenon a dual front. 
The aim of this paper is to explain this somewhat counterintuitive fact, and to put it into a more 
general perspective. Diffusion problems with several coexisting fronts have been discussed also 
elsewhere in the literature, see e.g., [PNJ], [GPP]. The mechanism described in those papers is, 
however, of a different nature from the one studied here, since the intermediate state is stable. 
Before starting the general discussion of dual fronts, we recall the model of [EBS]. 

In that paper, a one-dimensional ferromagnetic material is considered, which is described 
by the magnetization M(x) £ R 3 , with |M(a;)| = 1. Supposing that the material has an 
easy-plane orthogonal to the direction £, one obtains for the free energy 

W = i J dx (|M'| 2 + M 2 ) 

= ^ J dx (tf' 2 + <p' 2 sin 2 t? + cos 2 #) , 

where ' stands for derivative with respect to x, is the component of M in the direction £, 
and <p and i? are spherical coordinates with respect to polar direction £. It is assumed that the 
dynamics is driven by the Landau-Lifschitz equation [HC], i.e., 

M = - 7 (M x H) — \M x (M x H) , 

where H = —SW/SM, and that damping is sufficiently high such that the torque part 
— 7(M x H) can be neglected. Setting A = 1 yields the evolution equation 

# = #" + (!_ ) sin cost? , 
ip = ip" -\- 2(p' , &' cot i? . 



(1.2) 



(1.3) 
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In terms of the local wave number k = <p' , the second equation in (1.3) becomes 

k = k" + 2(&(lnsini?)')' . (1.4) 

We emphasize that the singularities at i? = mr, n £ Z are artificial, as can be seen by considering 
the complex variable A(i?, </?) = sini?e l¥ \ Equation (1.3) is equivalent to 

A = A" + ^-\Af + \Ar + ^^j A, (1.5) 

which is regular around A = 0. 

The following uniform stationary states of (1.3) will be considered: 

1) i?(a:) = 7r/2, <p(x) = k Q • x , 

2) $(x) = . 

Solutions of the first type are called by [EBS] spiral states, we call the second one the zero- 
solution. Linear stability analysis shows that spiral solutions are stable for \k \ < 1 and have 
one unstable direction pointing in the ^-direction when \k \ > 1. The solution $(x) = is 
unstable. 

If we choose an initial condition that connects a stable state to an unstable spiral state, we 
expect that a front will invade the unstable spiral. In [EBS] this phenomenon was studied in 
detail, for typical initial conditions connecting the stable state with k = to an unstable spiral 
state, i.e., for 

7T 

lim $(xit = 0) = — , lim <p (x,t = 0) = k Q > 1 , and lim <p (x,t = 0) = . 

x — >±oo 2 x — >oo " x — ► — oo 

It was observed that for k > \[2 the solution converges towards a leading front that connects 
the unstable spiral solution to the zero- solution, which in turn gives rise to a transition to the 
stable solution with <p'(x) = 0. The size of the region in which A is close to zero was seen 
to grow with time; see Fig. 1. This phenomenon of a growing intermediate region is the main 
object of this paper. 

Furthermore, a third front, which we call "ghost front," can be seen in the simulations. It 
is the line going through zeros of the amplitude \A\, which is visible in Fig. 2. Its origin will be 
explained in Chapter 5. 

The structure of this paper is as follows. We will show that the phenomenon described in 
the magnetic model above is a special case of a more general set of similar phenomena which 
can be observed in many parabolic models. After a discussion of fronts in Section 2, we will 
start by exhibiting this phenomenon in 1 -component systems in Section 3, then we describe a 
structurally stable situation with 2 components in Section 4, and finally describe the magnetic 
model in Section 5. Most of our results cannot be rigorously proven, but a few facts can be 
established rigorously. Some of them are collected in the Appendix. 
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Fig. 1: Several time-frames (t = t Q to t = t Q + 3) for the evolution of Eq. (1.3), with k Q = 2. One can clearly see the 
leading front moving with a speed c 12 and the trailing front moving with different speed c 23 < c 12 . 




Fig. 2: Several time-frames (t = t Q to t = < + 3) for the evolution of Eq.(1.3), with k Q = 2. One can clearly see 
the leading front moving with speed c 12 and the trailing front moving with speed c 23 . In the center, one recognizes the 
ghost front, moving with speed c ghos( . The vertical scale has been chosen in such a way to make the phase-slip points 
in the ghost front better visible (these are very small amplitudes). The three solid lines indicate the speeds. The middle 
line is fitted through the minima of the amplitudes. The boundaries of the "central valley" expand with the speeds c* l2 
and c 23 , since the overall amplitudes are going to zero with the same speeds as the trailing parts of the two fronts. 

2. Fronts and their speeds 

Before we can discuss the existence of dual fronts, we need a precise description of such notions 
as "critical speed" and "minimal speed" for fronts of parabolic equations. This discussion will 



-5- 



follow mostly the original definitions of Aronson and Weinberger [AW]. Let us consider for 
simplicity an equation of the form 

u = u" + f(u), (2.1) 

where dots stand for temporal and primes for spatial derivatives, and where u(x,t) is a real 
function. We assume that / £ C 1 , and that, on the interval [£,r], the function / satisfies 
f(x) > for x £ (£, r), and f(£) = f(r) = 0, see Fig. 3. 



f(u) V(u) 




> u 



Fig. 3: Typical shape of the force / and of the corresponding potential V. 

As is well known [AW], the existence of fronts (travelling waves) and their speed c is 
discussed by going to one (or several) moving frame(s) and studying solutions of the form 
u(x,t) = q(x — ci). Then q satisfies the "Hamiltonian" equations 

q" = ~cq' ~ V'(q) , V(q) = f(q) , (2.2) 

or, equivalently, 

\ =? ' (2-3) 
p = - cp - f(q) . 

One sees that c appears as a friction, and / as the derivative of a potential V in which the 
particle is moving, so that the zeros of / correspond to equilibria of this potential. We begin by 
discussing the existence of fronts, in the case of the potential V with critical points at I and r. 
(See also similar discussions in [CE], [AW], [PNJ].) Using the picture of a Hamiltonian flow 
with friction c, we observe that when c = 0, a particle released at r will "fall down" and traverse 
the critical point I ("overshooting"). 

Definition. The minimal speed associated with / is the smallest friction c for which the particle 
starting at the point r with zero velocity reaches the point I without overshooting. We call this 
speed c£ in . 
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Remark. In the case under discussion, the existence of c™ n is obvious. However, in the next 
section we will see an example where certain minimal speeds do not exist. 

Definition. The critical speed c* lv is the infimum of those frictions for which there is a initial 
velocity at r (which may be zero or non-zero) allowing a particle leaving r to the left to reach I 
without overshooting. 

Remark. The critical speed always exists. See Appendix A for a precise definition of c* lr . We 
also state there that c™ n = c* lr for the case considered in this section. 

An important issue in the theory of fronts is their "selection," i.e., the question which of 
several possible speeds will be chosen by an arbitrary initial condition. One has the following 
fundamental theorem about front selection that was proven in [AW]: it shows that the dynamics 
chooses the front with speed c* lr . Assume that / is as above. 

Lemma 2.1. [AW], [CE]. Assume that the initial condition u(x,0) takes values in [£,r], and 
reaches the limits 

lim u(x, 0) = I , 

x — >oo 

lim u(x, 0) = r , 

x — ► — oo 

for finite x. Then, one has 

lhn„(« + «M)={;; (2.4) 

Thus, if we move slower than the critical speed, we are "behind" the front and will see 
eventually the state r, and if we move faster, we see I. 



3. A one-component example 

We now discuss a one-component model which exhibits a dual front. This model is again given 
by equations of the form 

u = u"+f{u), (3.1) 

where we now suppose that / is a non-negative C 1 function, f(u) > in [0, 1] with three zeros 
at S 1 = 0, S 2 = |, and S 3 = 1; cf. Fig. 4. The choice of S 2 = \ is inessential. 

Coming back to the discussion of critical and minimal speeds, we can consider the three 
intervals (£, r) = (5 l5 5 2 ), (5 l5 5 3 ), and (5 2 , S 3 ). In the intervals "12" and "23", the discussion 
of the previous section applies without change: The minimal speeds exist, and are equal to the 
critical speeds, which we denote c\ 2 and c^ 3 , respectively. The critical speed c\ 3 for fronts 
connecting S 1 and S 3 exists by definition, but there is not necessarily an orbit with minimal 
friction connecting S 3 with 5 1 , since the potential may be such that any friction which lets the 
particle traverse S 2 must be so weak that the particle will necessarily overshoot at S 1 . In fact, 
this is the case for the potential of Fig. 4. On the other hand, and this is the main point of our 
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Fig. 4: Typical shape of the function / and of the potential V. 




choice of V, non-existence of c™ 3 m is equivalent to the inequality c* 12 > c 2i . As we will see, the 
strict inequality, c* 12 > c* 2i , implies the phenomenon of a dual front. 

We apply the Lemma 2.1 to the two subintervals "12" and "23". This will prove that under 
the conditions mentioned above, the size of the intermediate state grows linearly in time for all 
initial conditions that approach the limiting values sufficiently fast. We have a dual front: 

Proposition 3.1. Let f be such that c* 12 > c* 2i . Then, for any initial condition taking values in 
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[S 1 , S 3 ] and reaching the limits 

lim u(x, 0) = S 1 , 

x — >oo 

lim u(x, 0) = S 3 , 

£ — ► — OO 

for finite a;, an observer moving with speed c will see the limit 

if c* 3 > c , 
II c 12 > c > c 23 , 
if C > Cj 2 . 

Remark. The result above can be interpreted as follows: If we look in the laboratory frame at 
a position ct, with c £ (c 23 , c^ 2 ), we will find u very close to S 2 if is sufficiently large. Thus, 
the state S 2 is visible on an interval whose length is of order (c* 12 — c 23 )t. On the other hand, 
at x = ct, with c > c\ 2 , u will be close to S 1 . At x = ct with c < c 23 , we see the state S 3 . 

Proof. Case c* 12 > o c 23 : 

Consider the functions ip(x, 0) = min('u(a;, 0), 5 2 ) and ^(a;, 0) = max('u(a;, 0), S 2 ) as initial 
data for the Eq.(3.1). The Maximum Principle [F] tells us that 

ip(x,t) < u(x,t) < ip(x,t) for all (x,t) G R X R^" . 

Moreover, Lemma 2.1 tells us that lim^^ ip(x + ct, t) = S 2 for all c > c 23 and that 
lim^^ ip(x + ct, t) = S 2 for all c < c\ 2 , which completes the proof in this case. 

Case c > c\ 2 \ 

The proof for the case c > c^ 2 , is slightly more complicated. We consider a sequence of 
functions f s converging to /, as 8 { 0, such that f s (x) > for all x £ 5 3 ). We choose 
them such that f(x) + 8 > f g (x) > f(x) and f s — f has support in [S 2 — 8, S 2 + 8]. Thus, the 
corresponding potential V s has only critical points at S 1 and S 3 and Lemma 2.1 can be applied 
to the interval [S 1 , S 3 ]. The critical speed cl 3 (8) is equal to the minimal speed c™ 3 m (<5>). The 
proof for the case c > c\ 2 will follow from 

limc* 3 (£) = c* 2 . (3.2) 

o — 

To show Eq.(3.2), we first observe that f s (x) > f(x) implies cl 3 (8) > c\ 2 , by monotonicity, 
cf. also Appendix A. On the other hand, for every c > c\ 2 , there is a sufficiently small 
perturbation f s of / for which a particle released at S 3 with zero velocity will actually arrive 
at S 1 without overshooting. To see this, note that for any c > c\ 2 , there exists a /i(8) > for 
which a particle starting at q = S 2 — 8, > p > —fj>(8) will reach S 1 without overshooting. 
(This follows from the definition of c* 12 as the limit of such initial data.) By continuity, there 
is a /j, > such that fj,(8) > /j, for sufficiently small 8 > 0. On the other hand, if 8 > is 
chosen sufficiently small, the orbit in the potential V s starting at S 3 with zero speed will reach 



lim u(x + ct, t) 



^3 5 
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a point q = S 2 — 8, p = —u, with v > as small as we wish. As soon as v < /j, , this means 
that the particle will reach S 1 without overshooting, showing that cl 3 (8) < c. Since this is true 
for all c > c* 12 , the proof of this case is complete. 



Since the solution u satisfies u(x,t) G [S-^Sg] for all t > 0, we can change / for x < S x 
without affecting the solution u. We will do this in such a way as to make the critical damping 



from the left side less than c 23 . Define 



f(x) , ifx>S lt 
A • (x — 5 1 ) , if x < S 1 . 



We choose < A < (c2 3 /2) 2 . Next, we fix d satisfying 

2a/A < (Z < c* 3 . 
Consider now the "evolution" equation 

q ~ P ' (3.3) 
p' = -d-p- f x (q) . 

If we release a particle at S 3 with zero speed it will — due to the definition of c^ 3 and the choice 
of d < — overshoot at S 1 . Since the potential associated to f x is increasing (to 00) as 
x — >■ — 00, the particle will eventually return to S 1 from the left. This time, it will not cross 
S 1 again, because we have chosen d > 2^/\. Thus we have constructed a front for the force 
/ A , travelling with speed d which connects S 3 to S 1 . This front is not monotone, it overshoots 
exactly once; we have argued earlier that a monotone front does not exist in this case. Let us 
denote g X: d( x ) m e front obtained in this way — it is unique up to translation. By choosing a 
suitable translation we may and will assume that g x d (x) < u(x, 0). Applying the Maximum 
Principle to the problem u = u" + f\(u), we see that 

9x,d( x ~ td ) < «0M) > 
for all x and t. Thus, u(x + td' , t) > g x d (x + t ■ (d' — d)), for all d' , and hence 

lim u(x + td' , t) = S 3 , 

t — >oo 

when d' < d. Since this holds for every d satisfying 2 VA < d < c 23 , the assertion follows. We 
have covered all cases and the proof of Proposition 3.1 is complete. 

Discussion. We have seen that in one-component problems the situation is well-understood and 
that explicit conditions for the occurrence of dual fronts can be given. However, the reader 
should be aware that we are dealing here with a situation which is generically unstable with 
respect to changes in the potential. Namely, a generic change of / above will destroy the 
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critical point at S 2 so that it either disappears or splits into a maximum and a minimum. (This 
mechanism has recently been used by [PCGO] to characterize the critical speed.) Thus, the 
example is fully understood but not generically stable. In the next section, we consider a model 
where the order parameter has two components, and we will see that the analog of the saddle 
point S 2 remains structurally stable, so that the phenomenon becomes physically robust. On 
the other hand, rigorous proofs can only be given for part of the statements made, while some 
statements must be based on numerical experiments and a general analysis of the flow for the 
dynamical system. 



(4.1) 



4. A two-component model 

To be specific, we consider the following explicit model: 

dV 

a = (1 — 8)a" + — — = (1 — 8)a" + (1 + e)a — a — aab , 
dV 

b = (1 + 8)b" + — = (1 + 8)b" + (1 - e)b - b 3 - aa 2 b , 

where 

V(a, b) = il + ^a 2 + { ^—^-b> - la 4 - Vfe 2 - . (4.2) 
v ' ' 2 2 4 2 4 v ' 

Here a > — 1, < e < 1 and < 8 < 1 are free parameters that will be fixed later. In the 

discussion below, we will not use that the Eq.(4.1) derives from a potential. However, it will 

be intuitively slightly more appealing to have a potential V, so that we can think in terms of 

friction driven motion, as in the previous sections. 

This model has four constant stationary solutions: 



(a,b) 



(0,0) = 2\ , 

(o,yi~7) = T 2 , 
(vT+7,o) = r 3 , 

(a 4 ,6 4 ) = T 4 , 



where 

1 + a l — a 1 + a 1 — a 

The phenomenon of a dual front will be observed for a transition from T x to T 3 going through 
the intermediate state T 2 . For this to happen, the unstable manifold of the state T 2 with respect 
to constant perturbations should have exactly dimension 1 in R 2 , and the final state T 3 should 
be stable. This is true if and only if 

1 - £ 1 + £ . , 

< a < . 4.3 

1 + £ 1 - £ v ' 

Furthermore, we need conditions which establish inequalities between the various front speeds in 
suchaway that the phenomenon of a dual front can take place. Thus, we want c 23 < c* 13 < c* 12 . 
For this to happen, we impose 

8 > £ . (4.4) 
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Conjecture 4.1. Under the conditions (4.3) and (4.4) one has the inequalities 

C 23 < C 13 < C 12 • (4-5) 

Some arguments in support of (4.5). First, we can show that c™ n = 2y/i — ey/l + <5>. Indeed, 
if we consider only orbits from T 2 to T x with a = 0, then the 1-component analysis applies, and 
we get therefore (in the 2-component context) c™ n < 2y/l — ey/l + 8. This latter quantity is 

2^/T^~8~(dlV\ a=b=0 ) . Because the unstable manifold at T 2 has dimension 1, there can be 
no other optimizing orbit, and the result follows. Similarly, 

c™ n < 2y/lT^VT^8 < c™ n , (4.6) 

by our choice of e and 8. But now, we cannot show that the first inequality in Eq.(4.6) is an 
equality, because the dimension of the stable manifold at T 3 is 2. 

Consider next a front connecting T 2 to T 3 . Here, we can only give a local analysis near 
the point T 2 . The flow near T 2 leads to a monotone function a(x) and thus, we can describe the 
pair (a(x — ct), b(x — ct)) by the graph b(a) describing the front near T 2 = (0, y/1 — e). We 
have 6(0) = (1 - e) 1/2 . Substituting in Eq.(4.1), we get 

a = (1 - 5)o" + g b {a) , 

with an effective force 

g b (a) = (1 + e)a — a — aa(b(a)) . 

We reapply the discussion of Appendix A to this force, and we get the inequalities (at least 
locally near a = 0), 




It is this argument which we cannot really make global, because of insufficient control over b 
and its dependence on the speed c. On the other hand, numerical studies show that in fact a 
is monotone and the inequalities seem to hold globally. Moreover, the selection mechanism 
chooses c 23 = 2y/(l — 8)g' b (0), at least as long as the local force is concave. Thus, the 
situation seems to be exactly the same as the one encountered in the rigorously controlled 
example described in Lemma A. 2. 

In Fig. 6 below, we show the results of a numerical simulation. The dual front is clearly 
visible. On doing numerical experiments, one sees that the inequalities (4.7) hold, and that a 
selection mechanism takes place, which is similar to the one we rigorously established for the 
1-component systems. Here, we do not have a proof. But we can illustrate what is going on in 
a more intuitive picture, see Fig. 7. 
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a(x, t) 




^ x 



Fig. 6: Time evolution for the model Eq.(4.1), for the parameter values a = e = 0.5, 8 = 0.8. The time-frames 
shown are t = 5 to t = 50. Starting from an initial condition with very little 6-content, a 6-front is seen to emerge, 
which runs faster (with speed c* n ) than the a-front (speed c^). Note that the b state decays to with the speed Cj 3 , 
where the a-amplitude grows. 




Fig. 7: The potential V as a function of a and b for the same parameters as in Fig. 6. T 3 is a maximum, T x is a 
minimum, and T 2 is a saddle. Superimposed are orbits for times t = 5, t = 12, t = 15, < = 18 and t = 50. Note that, 
as time increases, the orbit comes closer and closer to a "circuit" T 3 -T 2 -T l . 



We are really considering a particle in the potential V. Since the "masses" 1 — 8 and 1 + 8 
are different we can imagine for fixed c the motion of a point particle in the potential V with 
two different frictions c/(l — 8) > c/(l + 8) in the a and b direction, respectively. The critical 
dampings then satisfy the inequalities (4.5). 
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Discussion. The new feature in the 2-component system is the existence of a monotone front 
connecting the state T x to T 3 , whereas in the 1 -component case, such fronts (connecting S 1 to 
S 3 ) did not occur. Nevertheless, the phenomenology of the 1 -component systems survives, as 
shown by our numerical experiments, and illustrated in Fig. 6. The intuitive reason for this can 
be seen in Fig. 7. We see that typical initial data pick up an orbit coming closer and closer to T 2 
as time advances, so that the transition through the intermediate state T 2 is "more advantageous" 
than a "direct transition" from T x to T 3 . This phenomenon is stable under small perturbations of 
the equations, but harder to prove, whereas the 1-component case was unstable, but accessible 
to rigorous control. 

5. Dual fronts in the magnetic model 

In this section, we come back to the model of the easy-plane ferromagnet described in Eqs.(1.3) 
to (1.5). The discussion will now follow very closely that of the 2-component model discussed 
above, with a third "ghost front" caused by topological effects. In this model we have again 
three stationary solutions which will be connected by the fronts: 

M 1 : $ = tt/2 , <p' = k , 
M 2 : ■& = , 

M 3 : ■& = tt/2 , <p' = , 
with k > y/2. One can define the front speeds as before and one has 



c 2 * 3 = 2 , c* 13 = 0, c\ 2 = 2^k 2 - 1 > 2 . (5.1) 

This can be seen by the following arguments: 

- The front connecting M 2 to M 3 corresponds to a solution with <p'(x) = , and therefore the 
1-component discussion of Appendix A applies. Since the function sin i? cos i? is concave, 
this means that the speed is given by the derivative at i? = 0, i.e., c 2i = 2. 

- The front connecting M 1 to M 3 with ■& = tt/2 has speed because (1.4) reduces in this 
case to k = which is the diffusion equation. 

- The speed of the front connecting M 1 to M 2 is discussed in much the same way as in 
the 2-component model. We assume that i? and <p' are monotone functions of x, and 
consider the function &(i?) = <p'(i?). Again, we have a family of concave force terms, 
(1 — &(i?) 2 ) sin i? cos i?, and therefore the derivative at M 1 determines the speed, leading 
to c* 12 = 2-yJkl — 1, as asserted. 

Applying the arguments of the preceding sections, we see the appearance of a dual front 
with an intermediate state which gets larger as time advances. This situation is fully confirmed 
by the numerical simulations of [EBS], where this phenomenon was discovered. One can see 
this clearly in Fig. 1 and Fig. 2, where a numerical simulation for this model is shown. Note 
that there is a third front appearing between the leading front travelling with speed c* 12 and the 
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trailing front travelling with speed c 2Z . We will call it ghost front for reasons which will become 
obvious shortly. 

The behavior of the solution near the intermediate state M 2 is most efficiently analyzed in 
terms of the equation for A = sini?e l¥ ', cf. Eq.(1.5). Since M 2 corresponds to A = 0, we may 
linearize Eq.(1.5), leading to the trivial equation 

A = A" + A . 

The asymptotics of the trailing front leads to a factor A 1 e~ x+2t in the laboratory frame, while 
the asymptotics of the trailing part of the leading front is 

A ^ e \x + (1 + X 2 )t _ 

The constant A can be determined from the leading part of the leading front, by observing [EBS] 
that the r.h.s. of Eq.(1.3) only depends on the derivative of <p and thus leads to a conserved 
quantity C, given by 

d^k + c* 12 k + 2Jfed e (lnsin#) = C, 

where, £ = x — c* 12 t. Thus, the leading part of the leading front determines its trailing part, and 
hence A. The explicit expression for A is given in [EBS]. 

Coming back to the problem of the ghost front, we want to determine the position of the 
zeros of the sum A 1 e~ x+2t + A 2 e Xx+ ^ 1+x as a function of time. Solving for the zeros leads 
to the equation 

(A + l)x + (1 + A 2 - 2)t = log(-A 1 /A 2 ) + 2iirn , 

with n £ Z, so that if we have one zero at (a; „,<„), then the others are — in this linear 
approximation — at 

(ImA) 2 + 1 - (Re A) 2 1 + Re A 

*• + 2 ™ |l + ApimA ' <° + 2 ™ \l + X\Hr*X 

Thus, the zeros seem to move with speed 

, (ImA) 2 

c , . = 1 - Re A + '— . 

ghost 1 + Re A 

This third speed is clearly visible in Fig. 2. The prediction for c ghost made in [EBS, Eq.(19)] 
introduced some spurious higher corrections in A but differs by less than 2% from the one made 
here. They both agree very well with the findings from the simulations. 
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Conclusion 

In this paper we have investigated the splitting of fronts propagating into linearly unstable states 
(dual fronts). Here, a single propagating front defining an interface between an unstable state 
and a stable final state is either unstable or may not exist. Instead, two fronts propagating with 
different velocities will emerge. They build up a linearly unstable intermediate state whose size 
grows linearly in time. In the frame moving with the leading (i.e., faster) front the intermediate 
state is only convectively unstable but it is absolutely unstable in all frames travelling slower 
than the trailing (i.e., slower) front. 

This scenario can be rigorously proved for a structurally unstable 1 -component model. 
Furthermore, we have good control for a generic 2-component reaction-diffusion model. Finally, 
we can explain in the same way this phenomenon in the magnetic model [EBS] where it was 
first discovered. It turns out that a third front can be observed in this model: The phase-slip 
front which is a ghost front resulting from the linear superposition of the tail of the leading 
front and the head of the trailing one. Such ghost fronts in the intermediate state are generally 
to be expected if the stable state and the unstable state are spatially periodic with different 
wavelengths. 

The main mechanism responsible for dual fronts can be understood in the following 
intuitive way: Suppose we disturb locally and weakly the unstable state. In the early evolution 
the dynamics is governed by the linearized equation of motion. Thus the perturbation can be 
understood as a linear superposition of pulses, each of them governed by a linear equation of 
motion which decouples from all others (e.g., pulses in a and b in the 2-component model). A 
necessary condition for front splitting is the fact that the pulse which expands fastest leads to an 
unstable state from which another perturbation will drive the system into a stable state. This is 
the meaning of the inequality c* 12 > c^ 3 . 

Acknowledgements. This work was partially supported by the Fonds National Suisse. 



Appendix A. Definition of the critical speed 

We have given an "intuitive" definition of c* lr in Section 2. Here, we present its mathematical 
counterpart, as given in [AW], for / as in Section 2. For monotone fronts, we can rewrite 
Eq.(2.3) as 

^ = _ C _M. (A . 1} 
dq p 

Let p cs (q) be the solution of Eq.(A.l) in the strip (q,p) £ [£, r] x ( — oo, 0] with initial condition 
Pc,s(^) = ~ £ - If Pc,s(<l) = for some q > £, then we define p C:S (q) = for all q > q. 
Since these curves cannot cross, it follows that they are monotone in e, that is, e < e implies 
Pc^il) > Pcs'il) for all 5 £ [£,r]. Since p c e (-) < 0, it follows that the limit 

exists and is a curve lying in [£, r] x ( — oo, 0]. Finally, we define 

c* r = inf{c > | p c (r) < 0} . (A.2) 
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In some special situations it is possible to determine the value of c* lr very easily from the 
nonlinearity /, namely when 

f( x ) 

sup — — - = / {£) , 

and / is positive on (-£, r). Then, 

A more detailed statement can be based on the following definition. For simplicity, we assume 
I = 0, r = 1. We also denote c* = c^. 

Definition. Define d~ = /'(0), d + = sup 0<a!<1 f(x)/x, and then c^ 1 = 2y/d±. 

Then, in their fundamental work [AW], Aronson and Weinberger prove the following 
results: 

Proposition A.l. Let f(0) = /(l) = andf(x) > forx e (0,1). Then 

i) For any c > c* there exists an orbit P c of Eq.(A.l) connecting (q,p) = (0, 0) to (q,p) = 
(1,0), which stays in the strip (q,p) £ [0, 1] x ( — oo, 0]. 

ii) No such solution exists for c < c* . 

Hi) If c* > 2y/ f'(0), then the connecting orbit P c coincides with p c , if and only if c = c* . 
iv) Ifc>c~ then p c has near (0, 0) the slope \ c , where 

K = -|(c+^ 2 "4/'(0)) , (A.3) 
and every other curve through (0, 0) has slope \ c with 

K = -|(^-^ 2 -4/'(o)) . 

v) (Selection) The critical speed c* satisfies c~ < c* < c + . If we take a localized positive 
perturbation of a front as initial condition, then the solution to Eq.(2.1) will eventually 
travel with speed c* . The shape of the selected front is given by the curve P c , . 

Proof. The proofs can be found in the following places in [AW]. 

i) Theorem 4.1 and Remark 1. 

ii) Lemma 4.3. 

iii) Theorem 4.1 and Remark 1. 

iv) Proposition 4.4. 

v) Proposition 4.2, the definition of c*, Eq.(A.2), and Theorem 5.1 and 5.2. 
Remarks. 

- Since i) and ii) say that a connecting orbit exists for c > c*, but not for c < c*, this means 
that c mm = c*, as asserted in Section 2. 

- If /'(l) 7^ then (0, 1) has a one-dimensional unstable manifold for the flow defined by 
Eq.(2.3) and the orbit p c mentioned in i) is unique. 
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Fig. 8: The three curves c , c + , and the front speed c* which is only different from c for fi > 2 although c + and 



- The result v) above has the following consequence: If / is concave, then c* = c = 



- One can improve the upper and lower bounds by considering polygonal approximations 
to / instead of the linear ones used above. 

A.l. Bounds on the selected speed 

This subsection illustrates the problem that the selected speed of a front is difficult to determine, 
and does not follow from any simple functionals of /, such as /'(0) ormax I>0 f(x)/x. This in 
turn is responsible for the difficulty of giving bounds on the selected speeds for problems with 
more than 1 component. 

The question we ask is whether it is possible that c~ = c*, even when c~ ^ c + . This 
question has been answered positively by Hadeler and Rothe [HR]. We give a slightly different 
account of the argument of [BBDKL] . 

Consider the special case of f^u) = u + (/j, — l)u 2 — /j,u 3 . Define in analogy with 
c 1 * 1 , for the function f . We have the 

Lemma A.l. For f as above, the critical speed c* satisfies c* = 2 for /j, £ [0,2], but c+ > c~ 
for /j, > 1 . 

Remark. The various c's are shown in Fig. 8. One has 



already differ for fi > 1. 



277 7 (0). 
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Proof. Choose the function v(£) = \ (l — tanh( 7 £)). Then, some algebra leads to 



v" + cv' + f^v) = |(l-tanh( 7 C)) 2 7 2 tanh( 7 0-c 7 +| + f(l-tanh( 7 0) • (A.4) 



(The signs are the same in both expressions.) Since we are interested in c > 0, only the plus 
sign is interesting for us. The calculation above means that we have found a positive front 
solution v = v , moving with speed C ' to the right. Therefore c* < C . But we also know 
from Proposition A.l, v) that c* > 2, since c~ = 2 for all /i. Furthermore, for /i < 1, the 
function / is concave in u £ [0, 1], and therefore c* = 2 for /j, < 1. Since we have found a 
function v for /i = 2 which produces a front with speed C ^ = 2, we find c\ < 2 and therefore, 
in fact c\ = 2. We shall complete the proof of Lemma A.2 by showing that c* is a monotone 
(increasing) function of /j,. Thus c* = 2 for all /j, £ [1,2], and we get the desired result. 

To prove the monotonicity, observe first that f^(x) is a monotone increasing function of 
li, for x £ [0, 1]. So, for all c, and for all /j, < 2, one has 

for all q g [0, 1]. Therefore it follows from Eq.(A.2) that c* < c^, if n < 2. 
Next we show 

Lemma A.3. For / as above, the critical speed c* satisfies c* = for /j, > 2. 

Proof. First we note that for c = C ^ the connecting orbit ^ is a trajectory p c defined by the 
above mentioned e-limit. This can be seen by using Proposition A.l, iv) and simply verifying 
that v has near (0,0) slope A c , see (A. 3). This result implies together with Proposition A.l, iii) 
thatc; = C M . 

Remark. The proof really shows that if there is a monotone connecting orbit with speed 
which coincides with the "fast" stable manifold of the point (0, 0), then we have c* = . 
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